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ABSTRACT 

This paper presents a robust controller design for second-order dynamic systems. The controller is 
model-independent and itself is a virtual second-order dynamic system. Conditions on actuator 
and sensor placements are identified for controller designs that guarantee overall closed-loop 
stability. The dynamic controller can be viewed as a virtual passive damping system that serves to 
stabilize the actual dynamic system. The control gains are interpreted as virtual mass, spring, and 
dashpot elements that play the same roles as actual physical elements in stability analysis. Position, 
velocity, and acceleration feedback are considered. Simple examples are provided to illustrate the 
controller design. From this illustration, the physical meaning of the controller design is apparent. 

INTRODUCTION 

Control theory for time-invariant linear systems which are described by first-order dynamic 
equations have been well established for decades. Control software tools today are also written in 
first-order forms. For applications, engineers can simply convert whatever models they have to 
the first-order forms and then use the existing tools to design the controllers. If the performance 
requirements are satisfied by the controllers, the design jobs are completed. If not, the design 
parameters are changed and the design procedure continues until a satisfactory design is found. For 
a small scale system, a few design iterations may be enough to complete a satisfactory design. 
However, for a large scale system such as the space station, the dynamic model usually involves a 
large number of degrees of freedom and is best described by second-order dynamic equations in 
terms of sparse structural matrices including mass and stiffness matrices. For second-order 
dynamic systems, transforming to first-order form not only increases the dimension of the 
problem, but also destroys the sparsity of the structural matrices, i.e. the mass and stiffness 


matrices for flexible structures. As a result, computational efficiency and physical insight are lost 
in the first-order form. Existing control analysis and design software may not be able to handle 
such a large system due to computational difficulties. For example, solving a 1000- by- 1000- 
dimension Riccati equation is considered numerically impossible using today's numerical 
techniques. There are basically two ways to address the controller design problems for a large- 
scale system. One way is to minimize the dimension of the system model by first preserving the 
second-order form and then performing model reduction. Laboratory experiments are required to 
verify the reduced model for robust controller designs. Recently controller designs using second- 
order system equations directly have gained attention in the literature as identified in Ref. [1]. 
Their computational advantages and physical features are also illustrated in Refs. [2] and [3]. 
Another way is to design a model-independent controller, which is insensitive to system 
uncertainties. The objective of this paper is to derive model-independent controllers for dynamic 
systems using second-order dynamic equations. 

When a mass-spring-dashpot is attached to any mechanical system, including flexible space 
structures, the damping of the system is almost always augmented regardless of the system size. 
The parameters of the mass-spring-dashpot are arbitrary, model-independent and thus insensitive 
to the system uncertainties. To satisfy the system performance requirements, the parameters are 
adjusted using the knowledge of the system model. The more the system is known, the better the 
parameters of the mass-spring-dashpot may be adjusted to meet the performance requirements. 
However, no matter what happens, the mass-spring-dashpot won't destabilize the system because 
it is an energy-dissipative device. The question arises as to if there are any feedback control 
designs using sensors and actuators which behave like the passive mass-spring-dashpot. This 
paper is motivated by this question and the answer is very encouraging. 

A novel approach for control of flexible structures is developed using a controller which can be 
described by a set of second-order dynamic equations. Under certain realistic (practical) 
conditions, this method provides a stable system with an infinite gain margin. For better 
understanding, two major steps are involved in developing the formulation of the method. First, 
consider only the direct output feedback for simplicity, implying the absence of dynamics in the 
feedback controller. Conditions are identified in terms of the number and type of sensors and their 
locations to make the system asymptotically stable with an infinite gain margin. Second, assume 
that the feedback controller contains a set of second-order dynamic equations. It is equivalent to 
visualize an imaginary flexible body, i.e. the feedback controller, which is linked side by side to 
the real flexible body. In other words, two sets of second-order dynamic equations are coupled to 
generate a closed-loop system. Design freedom increases when the dimension of the controller 



dynamic equations increases. Conditions are derived for the design of a stable closed-loop system 
having an infinite gain margin. The method takes advantage of the second-order form of 
equations (instead of transforming to a first-order form) which provides an easy way of discussing 
and obtaining the stability margin and results in a considerable computational efficiency for 
numerical simulations. Comparisons between the active feedback and the passive mass-spring- 
dashpot are given through several illustrative examples. 

DIRECT FEEDBACK 

In the analysis and design of dynamics and vibration control of flexible structures, two sets of 
linear, constant coefficient, ordinary differential equations are frequently used 

Mx + Dx +Kx = Bu (1) 

y = H a x + H v x + H d x (2) 

Here x is an n x 1 state vector, and M, D, and K are mass, damping and stiffness matrices, 
respectively, which generally are symmetric and sparse. The nxp influence matrix B describes the 
actuator force distributions for the p x 1 control force vector it. Typically, matrix M is positive 
definite whereas D and K are positive semi-definite. In the absence of rigid-body motion, K is 
positive definite. Equation (2) is a measurement equation having y as the m x 1 measurement 
vector, H a the m x n acceleration influence matrix, H v the m xn velocity influence matrix, and H d 
the m xn displacement influence matrix. Note that Eq. (1) can be solved for the acceleration in 
terms of the displacement, velocity and control force to obtain a new measurement equation in 
place of Eq. (2). However, physical insight is lost in this approach to controller design. 

The measurement equation, Eq. (2), may be used either directly or indirectly for a feedback 
controller design. Here we will use direct feedback. Let the input vector u be 

u = -Gy = -GH a x- GH v x- GH d x (3) 

where G is a gain matrix to be determined. Substituting Eq. (3) into Eq. (1) yields 


(M + BGH a )x + (D + BGH v )x + (K + BGH d )x = 0 


(4) 


For simplicity, consider the case where H a =H^= 0. Assume that the number of sensors m is 
larger than the number of actuators p. Let the actuators be located such that the row space 
generated by B T belongs to the row space generated by H v . In other words, the actuators are 
located in such a way that the control influence matrix B can be expressed by 

B T = C b H v (5) 

where C b is a p x m matrix which may be obtained by C b = B T H V T (H v Hj )T Assume that the 
gain matrix G is computed by 


G = LL r C b (6) 

where L is ap xp arbitrary matrix. Substituting Eq. (6) in Eq. (4) and noting the assumption that 
Hd = 0 leads to 


Mx + (D + BLL t B t )x + Kx = 0 (7) 

For the case where p < m < n, BLL T B T is positive semi-definite and thus D + BLUB T is at least 
positive semi-definite for a positive semi-definite matrix D. As a result, the closed-loop system, 
Eq. (7), is stable if D + BLL T B T is positive semi-definite, or asymptotically stable if D + BLUB T 
is positive definite. For the case where D is positive definite, D + BLLTB T is positive definite 
which yields an asymptotically stable closed-loop system. This leads to a conclusion that, for a 
structural system with some passive damping, an output velocity feedback scheme with non- 
collocated velocity sensors and actuators may make the closed-loop system asymptotically stable 
with an infinite gain margin since L in Eq. (6) for determination of the gain matrix G is an arbitrary 
matrix, as long as the actuators are properly located satisfying Eq. (6). Note that, for colocated 
sensors and actuators, B T = H v . Without velocity measurements, the system damping cannot be 
augmented from direct output feedback alone. However, if there are actuator dynamics involved, 
the system damping may be augmented by direct displacement or acceleration feedback. See Ref. 
[3]. 


CONTROLLER WITH SECOND-ORDER DYNAMICS 


Assume that the controller to be designed has a set of second-order dynamic equations and 
measurement equations similar to the system equations, Eqs. (1) and (2) 



M c Xc + D c x c + K c Xc B c 


( 8 ) 


y c — BIqqXc H vc Xc H dc %c (^) 

Note that this is a set of imaginary equations which do not represent any physical system. In fact, 
this set of equations basically serves as a filter to shift the phase of measurement signals. Here x c 
is the controller state vector of dimension n c , and M c , D c , and K c are thought of as the controller 
mass, damping, and stiffness matrices, respectively, which generally are symmetric and positive 
definite to make the controller asymptotically stable. The n c xm influence matrix B c describes the 
force distributions for the m x 1 input force vector u c . Equation (9) is the controller measurement 
equation having y c as the measurement vector of length p, H ac the p x n c acceleration influence 
matrix, H vc the p xn c velocity influence matrix and H dc the p x n c displacement influence matrix. 
Again, all the quantities, u c , y c > and n c are imaginary and thus arbitrary which means that M c , D c , 
K c , H dc and H vc are the design parameters for the controller. 

Let the input vectors u, u c in Eq. (1) and in Eq. (8) be 

u=y c = H ac x c + H vc x c + H dc x c ( 10 ) 

u c =y = H a x+H v x+H d x (11) 

Substituting Eq. (10) into Eq. (1) and Eq. (11) into Eq. (8) yields 

Mi'xt + D t x t + KtXt — 0 ( 12 ) 


where 


M t = 


M -BH ac 
-B c H a M t . 


,D t = 


D -BH VC 

~B C Hy Dc 


r = K -BH dc 
‘ -B c H d K c 



If the design parameters, M c , D c , K c , H dc and H vc are chosen such that M t , Dt and AT, are positive 
definite, the closed-loop system, Eq. (12), becomes asymptotically stable. 
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DISPLACEMENT FEEDBACK 


For better understanding of the advantage of the controller having second-order dynamic 
equations, consider a special case where H a = H ac = H v = H vc = 0. To make K t symmetric, it is 
required that 

BH dc = H T d Bj (13) 


or 


[b -Hi) 


H dc 

Bj 


= 0 


( 14 ) 


For the case where the sum of the number of actuators, p, and the number of sensors, m, is less 
than the number of states, n, the left-most matrix of Eq. (14) is a tall matrix. Unless B is in the 
space spanned by H d T or vice versa, there does not exist any solutions for B c , H dc in Eq. (13). 
Assume that the number of sensors, m, is larger than the number of actuators, p. Let the actuators 
be located such that the row space generated by B T belongs to the row space generated by H d , i.e. 
the actuators are located in such a way that the control influence matrix B can be expressed by 

B T = Q b H d (15) 

where Q b h ap x m matrix which may be obtained by Q b = B T H d T (H d H d T )- J . Substituting Eq. 
(15) into Eq. (13) yields 


HjQ T b H dc = HjBj (16) 

Since Hj is a tall matrix for m< n, the only possible solution is 

Q T b H dc =Bj (17) 

For any given matrix H dc , this equation produces a B C T which makes the matrix K t symmetric, 
i.e.. 


K t = 


K -BH dc 
_-hIb t K c . 


or 


Kt - 


K -HjBj 
-B c H d K c 


(18) 


The next question is how to choose a matrix H dc which makes the closed-loop stiffness matrix K t 
positive definite. The matrix K t is positive definite, generally written as K, > 0, if and only if 


xj K t x t >0 (19) 

for any real vector x, except the null vector. Substituting the definition of K, and x t from Eq. (12) 
in Eq. (19) yields 


xjK, x t = x T (K - HjBjB c H d )x + (B c H d x - x c ) T {B c H d x - x c ) + xJ{K c - I)x c 
= X T{K - B HdcHfc B T )x + {Hi B T x - x c ) T (Hj c B t x - x c ) + xJ{K c - T)x c 


( 20 ) 


This equation is greater than zero if B c and K c are chosen such that K - BH dc H dc T B T and K c - 1 are 
positive definite. Note that this is a sufficient condition but not a necessary condition. To make 
Eq. (19) hold, K must be a positive definite matrix, i.e. K > 0, and B c must be chosen such that 
K - BH dc H dc T B T > 0. It implies that this controller may not be able to control rigid body motion 
since K in this case is only a positive semi-definite matrix, K ^ 0. To release the constraint 
condition, K - BH dc H dc T B T > 0, K must be increased by at least BH dc H dc T B T . In other words, the 
system must be stiffened which can be achieved by adding displacement feedback. 

Let the input force be 


u.=y c - Gy = H dc x c - GH d x (21) 

where G is a gain matrix to be determined. Note that the velocity feedback is not considered here. 
Substituting Eq. (21) into the system equation, Eq. (1), the closed-loop stiffness matrix, Eq. (18), 
becomes 


K + BGH d -HjBj 
•B c H d K c 


( 22 ) 


If G is chosen such that 


G - H&Bc 


(23) 


which from Eq. (13) results in 


BGH d = BHfc B c H d = HjBjB c H d 


The closed-loop stiffness matrix, Eq. (22), thus becomes 


K t = 


K + BH dc Hj c B T -BH dc 
-hW K c 


or K,= 


K + HjBjB c H d -HjBj 
-B c H d K c 


( 24 ) 


which changes Eq. (20) to be 


xjK t x t = x T Kx + (B c H d x - x c f(B c H d x - x c ) + xJ{K c - T)x c 

r (25) 

= x T Kx + (//£fl T x - x c ) {Hjfi T x - x c ) + xJ{K c - f)x c 

Since K c is a design parameter, the closed-loop system becomes stable as long as K c is chosen 
larger than /, i.e. x c T (K c - I)x c > 0 for any arbitrary vector x c . An obvious choice is K e = / 
where / is an identity matrix of dimension n c . However, this is not the best choice which will be 
discussed later. To this end, it is shown that a stable closed-loop system can be designed using a 
feedback controller with second order dynamic equations. The controller has an infinite gain 
margin in the sense that the matrices M c , D c and K c , which may be considered as the gain matrices 
for the controller state vector x c and its derivatives, can be as large as desired without destabilizing 
the system as long as they are positive definite and K c is larger than I. 

A little modification of the above design produces a better design which has physical meaning. 
Indeed, let 


B c = K C B C or B c = Kc l B c (26) 

where K c is assumed to be positive definite so that the solution for B c exists for any given B c . In 
addition, let the gain matrix G in Eq. (23) be slightly modified as follows 

G = Hdc B c (27) 


which, with the aid of Eq. (13), results in 



BGH d = BH dc B c H d = hJbJk c B c H d 


The closed-loop stiffness matrix in this case (see Eq. (24)) thus becomes 


K,= 


K + HjBj_K c B c H d 

-K c B c H d 


-H d B‘K c 

K c 


( 28 ) 


which in turn changes Eq. (25) to be 


xjK,x t = x T Kx + (B c H d x - x c ) t K c (B c H d x - x c ) 
= x T Kx + (//J c B r x - x c ) T K c (Hj^x - x c ) 


(29) 


This equation is obviously positive if K is at least positive semi-definite, i.e. K > 0. Does this 
design have any physical meaning? The answer is positive. Consider the special case where the 
controller is as large as the system in the sense that the number of system states n is identical to the 
number of controller states n c . Furthermore assume that all the states are directly measurable, H d 
= I, and there are n actuators collocated with the sensors, B = /. In this case, Qb = / (Eq. (15)), 
B c = H dc = K c (Eq. (17)) for B c = I, and G = K C (Eq. (27)), which yields from Eq. (28) 


K t = 


K+K c 

-K c 


-K c 1 

K c J 


(30) 


For a single degree of freedom (n c = n = 1 ), K t represents the stiffness matrix for two springs 
connected in series with spring constants K and K c . 


PHYSICAL INTERPRETATION 

For better understanding of the nature of the dynamic control designs developed here, they are now 
interpreted in physical terms. In this section, three illustrative examples will be shown, starting 
with a simple spring-mass system. 

EXAMPLE I: A simple spring-mass system with a single-degree-of-freedom controller 
Consider a single-degree-of-freedom spring-mass system, n c =n = 1, with displacement 
measurement of the system mass. The second-order controller for this case reduces to a virtual 
spring-mass-dashpot system connected in series with the system mass as shown in the following 
sketch 


SYSTEM CONTROLLER 

(A Real Spring-Mass System) (A Virtual Spring-Mass-Dashpot System) 
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Fig. 1 A simple spring-mass system with a single-degree-of-freedom dynamic controller 


Let the position of masses m and m c be measured form their equilibrium states. The equations of 
motion for the above system can be derived by applying a force to m and m c . The force applied to 
the system mass m in this case is the force F transmitted through the spring k c . This is precisely 
the control force applied to the system as given in Eq. (21) with Hj c = k c , Hj = 1, and G = k c . 
Thus the second-order control law is simply 


u = F = k c (x c -x) 


(31) 


where x c is computed from 


m c x c + d c x c + k c x c = k c x 


The equation of motion which describes the closed-loop behavior of the above system is simply 


m 0 

X 

. 0 m c . 

x c . 


0 0 
0 d c 




(32) 


The above equation verifies Eq. (12) with H vc = H v = H ac = H a = 0 (i.e. no velocity and 
acceleration measurements), and K t given in Eq. (30). The above set of equations is always stable 
for any m, k, m c , k c , and d c , and is asymptotically stable for any d c > 0. We now consider 
various special cases. 
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Case 1: For the controller without damping, d c = 0, the system reduces to two spring-masses 
connected in series. If k c is small, the control force given in Eq. (31) is small, thus the controller 
exerts little influence on the system. Mathematically, Eq. (32) becomes a set of two uncoupled 



equations of x and x c , and obviously little change in the response of the controlled system is 
expected from this controller. If, however, k c is large, (i.e. the virtual spring is stiff) the relative 
displacement between the two masses is small. Hence in the limit the two masses move together 
like a single mass m + m Ct and the natural frequency of the system is approaching 

< 33 > 

. As a result, for large k c , changing the design variable m c will affect the natural frequency of the 
closed-loop system according to Eq. (33) above. 

Case 2: For d c > 0, the system is always asymptotically stable (unless k c = 0, which as discussed 
before means no control). The energy flows from m to m c and is dissipated by the damper. Again, 
for large k c , the system can be approximated as 

(m + m c )x + d c x + kx = 0 (34) 


Introduce the notation 


m + m c 


- = 2 Co)„ 


m + m c 


- a>l 


Thus, 


r _ J dc 

2 1k(m + m c ) 


(35) 


The design variables in this case are d c and m c . Various choices of d c and m c will result in C > 1* 
C < 1, or £ = 1, which corresponds to the cases the closed-loop system is over-damped, under- 
damped, or critically damped, respectively. 


Case 3: For general values of k c , d c , and m c , the design can be thought of as a virtual vibration 
absorber. Let the system be excited by some unknown force FeM 1 , and the displacement of the 
mass m be denoted by x = Xe^ at + The typical objective of a vibration absorber design is to 
determine the values of k c , d Cl and m c such that the ratio 


7 = 


FeiW 
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is minimized over an interested range of excitation frequency coy. 

EXAMPLE 2: A two-degree-of-freedom system with a single-degree-of-freedom dynamic controller 
Next, consider a two-degree of freedom spring-mass system with displacement measurements of 
the masses mj and m 2 from their equilibrium positions, n = 2. First consider the case where the 
the controller has only one state, n c — 1. The second order controller in this case is simply 
equivalent to a virtual spring-mass-dashpot system connected in series with the two system masses 
as shown below 


SYSTEM CONTROLLER 

(A Real Spring Mass System) (A Virtual Spring-Mass-Dashpot System) 



I J_ I 

Fig. 2 A two-degree of freedom system with a single-degree-of-freedom dynamic controller 


It can be easily shown that the control force applied to the system is simply 


.r^iU 0 


Fil IA te • xi) 


where Fj denotes the force applied to mj,j =1,2; and x c is given by 


m c x c + d c x c + k c x c = k c x 2 


Furthermore, the closed-loop behavior of the above system is governed by 


mi 0 Opt] [0 0 Olpil [*i+*2 -*2 0 Ipti] 

0 m 2 0 X 2 + 0 0 0 x 2 + -k 2 k 2 +k c -k c x 2 =0 (37) 

.0 0 m c JL.x c J l 0 0 z/ c J LJC 3 J L 0 k c JL^c- 
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which again verifies Eq. (12) with H vc = H v = H ac = H a = 0 and with K t given in Eq. (30). Note 
that the above scheme requires only displacement measurement of the mass m 2 - 


EXAMPLE 3: A two-degree-of-freedom system with a two-degree-of-freedom dynamic controller 
Consider the two-degree-of-freedom system above again with displacement measurements only, 
but now displacement measurement of the mass m\ is also to be used in the controller design. The 
second order controller design in this case is simply 


u = 


\Fi] 

Fi\ 


k Ci 0 j ' 
0 k C2 


Xc x - X\ ' 
X Cl - X2 . 


(38) 


where 


m c , 0 


0 m 


C2 


ic 2 j 


d Cx 0 
0 d, 


C 2 


'■Cl 


k Ci 0 

0 t 


Cl J 


'■Cl 

L^c 2 


_’k Cl oipn 
■[ o k Cl \[x 2 l 


(39) 


The closed-loop system is equivalent to a mass-spring-dashpot system shown below 


SYSTEM CONTROLLER 

(A Real Spring-Mass System) (A Virtual Spring-Mass-Dashpot System) 


r 



Fig. 3 A two-degree of freedom system with a two-degree-of-freedom dynamic controller 


whose behavior is governed by 
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( 40 ) 


m\ 0 0 0 

*1 


0 0 0 0 

*i 


k\ +kj +k Cl -ki -k c , 0 



0 m 2 0 0 

*2 

4- 

0 0 0 0 

*2 

+ 

-kj ki +k c 2 0 -k C2 


X2 

0 0 m C] 0 



0 0 d Cl 0 

* Ci 


l 

o 

o 


X Cl 

0 0 0 m Cl 

*C 2 _ 


. 0 0 0 d Cl 

1 

K* 

l 


1 

0 

1 

to 

o 

r 


IXC2] 


If velocity measurements are available, say at the system mass nt\, then a dashpot element may be 
added in between m\ and m c \ for example. It should be noted, however, that the controller 
masses, springs, dashpots are in fact virtual elements with physical interpretations as such. For 
ground based systems, they may represent actual physical elements attached to the ground. But for 
space based systems, they are simply controller gains in the control algorithm. 


ACCELERATION FEEDBACK 


The above controller can be extended to acceleration feedback as well. Consider the system given 
in Eq. (1), but now the measurement vector y in Eq. (2) has only acceleration measurements, i.e. 
H v = Hj = H vc = Hd C =0 in Eq. (12) 


M t x, + D t x t + K, x t = 0 


where 


M t - 


M -BH ac 
-B c H a M c . 


, D,= 


D 0 

L0 A:J 


K,= 


K 0 
L0 K c \ 


To make M t symmetric, it is required that BH ac = H a T B c T as discussed in Eqs. (13)-(17). All the 
discussions regarding the positive-definiteness of K t from Eqs. (18)-(20) also apply to M t . 


Additional coupling in the closed-loop mass matrix M t can be achieved by letting the input u in Eq. 
(12) include direct acceleration feedback, i.e., 


u = y c - G a y = H ac Xc - G a y 


(41) 


which makes M t become 


M,= 


M + BGH a 
. -B c H a 


- BHac ' 

M e . 


(42) 
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As before, M t can be made symmetric and positive definite by proper choices of H ac , B c , and G a . 


Let 


B c = M C B C or B c = Mc X B c (43) 

where M c is positive definite so that the solution to B c exists for any given B c . Let G be chosen 
such that 


G = H ac B c = H ac Mc X B c (44) 

which, with the aid of the equality, BHqc = HjBj, results in 

BGH a = BH ac B c H a = HjBjBc H a = HJbJm c B c H a (45) 

The closed-loop mass matrix in this case becomes 

M Jm + hJb[m c b c Ha -hJbJm c 1 (46) 

' [ -M c B c H a Me J 

This is a positive definite matrix as discussed in Eq. (29) for K t , regardless of the value of M as 
long as M is positive definite. The closed-loop in this case becomes 

M + hJb[m c B c H a -hJbJm c 

-M c B c H a M c 

Figure 4 shows a block diagram of the closed-loop system with acceleration feedback. In order to 
have a different flavor for the readers, all the quantities are expressed in frequency domain to 
quantify the multi- variable stability margins and performance of such systems. Let Gs(s) = H a 
[Ms 2 + Ds + K] X B be the system transfer function, G c (s) = Hac[M c s 2 + D c s + K C \ X B C the 
controller transfer function, and G a = H ac M c l B c the direct acceleration feedback gain. The 
acceleration measurement y(s) caused by the application of an external force r(s) can be expressed 
by 


x 

Xc. 


D 0 
0 D c 


x 

Xc 


o j£]L?J =0 <47) 


y(s ) = s 2 G s (s)[r(s)+(s 2 G c (s) - G fl )y(s)] 


(48) 
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or 


y(s) = [I - s 2 G s (s) (. s 2 G c {s ) - G*)]' 1 j 2 G s (s)r(s) 

The closed- loop transfer function from r(s) to y(s) is 

G(s) = [I - s 2 G s (s ) (s 2 G c (s) - G a )]-^ 2 G s (s) (49) 



Fig. 4 Block diagram of the closed- loop system with acceleration feedback 
It is interesting to note that 

s 2 G c (s) - G a = -H ac (/ Mcs 2 + DcS + ^ + K c ) M C X B C (50) 

Figure 5 is equivalent to Fig. 4 for the closed-loop system with acceleration feedback. All the 
quantities M c , D c , and K c are design parameters which are model independent but they must be 
positive definite. The quantities Hoc and B c are related by H a and B such that BHac = HffBj. This 
system is always stable regardless of how much uncertainties occur in the system matrices M, D, 
and K. 
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Fig. 5 Simplified block diagram of the closed-loop system with acceleration feedback 

Let Hoc be H ac = HacM c , and recall that B c = M c B c . Equation (50) becomes 

s 2 G c (s) - G a = -Hoc (s 2 + MjDcS + M c 1 K^\DcS + K c ) B c (51) 

For the case where M c is sufficiently large such that M C ~ 2 D C and M C ^K C may be neglected, the 
above equation can be approximated by 

s 2 G c {s) - G a = -Hoc {DcS + K c ) BcS- 2 (52) 

Figure 5 can then be reduced to yield Fig. 6. For the case K c = 0, the controller becomes an 
integrator of the acceleration measurement. If H ac is chosen to be B C T , then 

s 2 G c (s) - G a ~ - BjD c B c s 1 

which is equivalent to a direct velocity feedback to the system. 



Fig. 6 Reduced block diagram of the closed-loop system for a sufficiently large Me 
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Although it seems logical to choose a large mass matrix M c for the controller, measurement bias 
and noises may prevent such a choice in practice, because integrating a bias is obviously not 
desirable in a control loop. 

The procedure for deriving the second-order controller with acceleration feedback is identical to 
that for displacement feedback. Mathematically, both controllers are identical in the sense that the 
closed-loop mass matrix M t for acceleration feedback can be obtained by replacing K by M in the 
closed-loop stiffness matrix K t for displacement feedback, and subscript d by a. In other words, 
both displacement and acceleration feedback are conceptually dual. However, significant 
differences between both controllers appear when they are implemented either actively, or 
passively, which will be shown in the following example. 

EXAMPLE 4: A single-degree-of-freedom system with acceleration feedback 
Consider a single degree-of-freedom spring-mass system with acceleration measurement of the 
system mass. The second-order controller for this case reduces to a virtual spring-mass-dashpot 
connected in series with the system mass as shown in Fig. 7 


SYSTEM 


A 

^ — MM — | 
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-MM — I 

— E — 


m, 


Fig. 7 A simple spring-mass system with acceleration feedback 

Note that the vector x c here means the relative position of m c to the position of m. In this case, H a 
and B c in Eq. (47) are chosen to be B c = -H a = 1. The second-order control law is 

u = -m c (x c + ic) 


where x c is computed from 
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The closed-loop system can then be rewritten as 


' m + m c m c ' 

X 

+ 

' 0 

o' 

X 

+ 

k 

0 

X ' 

m c m . 

x c . 


. 0 

d c . 

X c . 


. 0 

k c . 

x c 


= 0 


The transfer function , Eq. (49), becomes 


G(s) = 


(ms 2 + k) + m c s^dc s + k c ) (m c s 2 + d c s + k c ) 1 


For large m c , G(s) reduces to 


ms 2 + dcS + (k + k c ) 

The system is clearly asymptotically stable. The numerator s 2 appears due to the acceleration 
feedback. 

Comparison of Fig. 1 and Fig. 7 reveals the difference between the acceleration and displacement 
feedback controllers. The controller for acceleration feedback does not have a virtual ground 
attached to the control mass and thus cannot control the rigid body motion. 

CONCLUSIONS 

This paper formulates a robust second-order dynamic stabilization controller design for second- 
order dynamic systems. The design is passive in the sense that it contains mechanisms that serve 
only to transfer and dissipate energy of the system. The controller interacts with the physical 
system only through spring, mass, and dashpot elements, and therefore, it can be implemented 
actively or passively. In other words, stabilization can be accomplished either by a controller with 
gains interpreted as virtual mass, spring, and dashpot elements, or by actual physical masses, 
springs, and dashpots connected to the system. 

The passive design means that the controller does not destabilize the system. As far as stability is 
concerned, the controller is model independent, and this is a robust design. Specifically, overall 
closed-loop stability is guaranteed independently of the system structural uncertainty and variations 


in the structural parameters. It should be emphasized that this is a robustness result with respect to 
structural uncertainty in the absence of measurement uncertainty and other contributing factors. 

However, control performance, unlike stability robustness, is dependent on the system 
characteristics. Knowledge of the system model can always help improve a controller design. In 
this method, the controller order and/or controller gains can be adjusted to meet the desired 
performance. Physical interpretation of the controller gains as virtual masses, springs, and 
dashpots provides convenient rules of thumb as to how they should be adjusted to meet a certain 
desired performance objective. 

Finally, the controller has been formulated from the continuous-time setting. Actual implementation 
of the controller, however, most likely requires usage of a digital computer. In future work, effects 
of sampling and time delays will be addressed. Other practical issues that can also affect the control 
performance such as measurement noises, actuator and sensor saturation limits will be 
investigated. It should also be interesting to examine the difference and relation between this 
approach and others such as full state feedback with a second-order state estimator. 
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